Effective Mass Anomalies in Strained Si Thin Films and Crystals 
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Effective mass anomalies due to the geometrical effects are investigated in silicon nanostructures 
using first-principles calculations for the first time. In {111} and {110} biaxially strained Si, it is 
found that longitudinal effective mass is extraordinarily enhanced for both thin films and crystals. 
This mass enhancement is caused by the change of the band structure with double minima into that 
with a single minimum due to strain and confinement. At the transition point, it is analytically 
shown that the effective mass diverges. The dependences of the confinement thickness on the 
anomalies are qualitatively explained by an extension of the effective mass approximation. 



In this decade, various nanostructures have been fab- 
ricated and measured. In such small structures, the ge- 
ometrical effect due to strain, confinement, and crystal 
directions significantly changes the electronic structures 
compared with that of the bulk. Intentional use of the 
geometrical effect has been suggested for the improve- 
ment of the device performance PJ. From the techno- 
logical viewpoint, the function of sub-nm-body metal 
oxide field-effect transistors (MOSFETs) using silicon 
on insulator [2j has been demonstrated, which consist 
of a silicon transport region sandwiched by insulating 
layers [3]. The nanostructure most sensitive to the ge- 
ometrical effect may be an ultrathin body double-gate 
(UTBD) MOSFET0, % with a body thickness of less 
than lOnm; it has a three dimensional structure with the 
channel region sandwiched by two gate insulators and is 
a promising sub- 10 nm device. Recently, in terms of elec- 
tron mobility for the UTBD MOSFET, the advantage of 
the (110) direction due to the valley structure has been 
reported in comparison with the (001) direction @. It 
is very interesting to investigate the change of the elec- 
tronic structures including transport parameters such as 
effective mass from the bulk to the thin films of some 
atomic layers. However, there are few studies on the 
intrinsic properties of semiconductor films in the inter- 
mediate region from the bulk to the sub-nm body using 
first-principles calculations Q- 

In this Letter, the effective mass is systematically in- 
vestigated for the biaxially strained Si thin films confined 
in the (001), (111), and (110) directions and for strained 
Si crystals using density functional calculations, where 
the strain is uniformly applied to the plane perpendicu- 



lar to the confinement direction. In the {111} and {110} 
strained Si, it is found that in some region, longitudinal 
effective mass is extraordinarily enhanced for the thin 
films and the crystals. Preliminary results are reported 
in The confinement and the symmetry breaking due 
to the strain cause these anomalies, where the band struc- 
tures change from double minima to a single minimum. 
In such a case, it is analytically shown that the effective 
mass diverges. As far as the author knows, this is the 
first report on the qualitative change of the band struc- 
ture and the divergence of the effective mass in Si caused 
by confinement and strain The dependences of the 
confinement thickness on the anomalies are qualitatively 
explained by extending the effective mass approximation 
to the whole band near the Brillouin zone boundary. The 
experimentally observable regions of the strain and thick- 
ness are also obtained for the (110) confinement. 

Band calculation is based on the density functional 
theory 1 id . Illj . The local density approximation (LDA) 



PW92[12j is adopted for the exchange correlation func- 
tional. To describe ion-electron interaction, norm- 
conserving [l3| and ultrasoft [lij pseudopotentials are 
adopted for Si and CI, and for H, respectively. The cal- 
culational condition is carefully checked [15j. The cutoff 
energy is 20.25 Rydberg. The meshes of the sampled k- 
points are 8 x 8, 8 x 8, and 4 x 8 for the (001), (111) and 
(110) slabs, respectively. The calculated longitudinal and 
transverse masses for the bulk without strain are 0.950 
(0.9163) and 0.189 (0.1905), respectively. The figures in 
the parentheses are experimental ones[16]. The unit of 
effective mass is the bare electron mass. The effective 
mass is calculated from a fitting of eigenvalues to the 6- 
th-order polynomial. The calculation is performed using 
Tokyo Ab initio Program Packages (TAPP) code [171 ]. 

As a model of Si film sandwiched by insulating layers, 
hydrogen-terminated Si (001), (111) and (110)-(1 x 1) 
slab models are used for the [001], [111] and [110] con- 
finements, respectively. The plane perpendicular to the 
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FIG. 1: Schematic illustrations of the electron pockets for the 
bulk and the confined thin films without strain. For the bulk, 
six equivalent ellipsoidal electron pockets are located along 
the r — X lines. Confinement projects these electron pockets 
onto the perpendicular plane: for (001), one isotropic valley 
(r~ valley) and four elliptic valleys (A- valleys); for (110) ,two 
kinds of inequivalent valleys (A- and A'- valleys), where the 
two and four pockets in the bulk correspond to the A- and 
the A'-valleys, respectively; for (111), six equivalent elliptic 
pockets. 



confinement direction is uniformly strained from -4 to 
+4% [20] and the lattice constant along the confinement 
direction is relaxed [18[- The slab models with hydro- 
gen and chlorine terminated for each side are addition- 
ally used to evaluate the surface effect. The Si atoms in 
the model consist of the crystal structure unit, which is 
optimized for the given strain. 

Before proceeding to the results, the electron pockets 
for the confined thin films are briefly introduced. For the 
non-strained case, Fig. [1] shows schematic illustrations of 
the electron pockets of the bulk and the slabs confined 
in the (001), (110) and (111) directions. For the (001)- 
and (llO)-confined slabs, there are two kinds of inequiv- 
alent valleys. The strain and the confinement thickness 
determine which valley is the bottom of the conduction 
band. Among them, the F- and the A-valleys are mainly 
focused. 

Fig. [2] shows the longitudinal effective mass at the A- 
valley in the [110] confinement as a function of the strain 
and the slab thickness. The strain is represented by the 
ratio of the lattice constant under strain to that for the 
strain-free bulk. The unit of thickness is nm. This figure 
clearly reveals that along lines, there are two regions with 
extraordinarily large masses. The confinement causes the 
contour lines to curve toward the direction of decreaseing 
strain. This tendency becomes stronger as the thickness 
decreases. For the [111] confinement, the longitudinal 
mass distribution is shown in Fig. 03 On the other hand, 
in the [001] -confined slabs, the behavior of the effective 
mass at the T-valley is very smooth and the variation 
range of the mass is narrow (0.18 ~ 0.20) in Fig. [H In 
particular, the effect of confinement is negligibly small. 




FIG. 2: Longitudinal effective mass at the A- valley (Fig. [I]) in 
the (110) -confined slab as a function of the strain and thick- 
ness. The unit of effective mass is the bare electron mass. The 
strain is represented by the ratio of the strained lattice con- 
stant to the bulk value without strain. The strain more/less 
than 1 corresponds to the tensile/compressive strain applied 
biaxially to the plane perpendicular to the confinement direc- 
tion. For the effective mass, values more than 50 are truncated 
to avoid spurious interpolation. 




FIG. 3: Longitudinal effective mass in the [111] -confined slab 
as a function of the strain and thickness. The configuration 
of the graph is the same as in Fig. [2] 

Note that, for the [111] and [110] confinements, a huge 
mass enhancement is observed in the strain-free region as 
well as the strained region. To evaluate the surface ter- 
mination effect, the slab models with the H and CI atoms 
for each surface side are investigated and it is found that 
these models also yield qualitatively the same results as 
above. 

It is reported that the behavior of effective mass for the 
slab with a thickness more than 10 nm is very similar to 
that for the bulkQ. Fig. [5] shows a comparison between 
the longitudinal effective mass at the A- valley for a slab 
with 65 Si atomic layer (~12 nm) (closed square) and 
for the bulk (open square) in the (110) strained case. As 
expected, both behaviors are very similar and reveal a 
mass divergence at the strains of 0.978 and 1.022. At 
these strains, the band structure qualitatively changes 
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FIG. 4: Effective mass at the T-valley (Fig. [Q) in the [001]- 
confined slab as a function of the strain and thickness. The 
graph configuration is the same as in Fig. [2] 



from double minima to an single minimum. In Fig. 
the k-position of the minima is plotted for the slab ( "+" ) 
and the bulk ("x"), where the position is represented by 
the fractional coordinates and the value 1 corresponds 
to the Brillouin zone (BZ) boundary. Therefore, the 
number of electron pockets for the strongly strained bulk 
changes from six to five. The same analysis is valid for the 
(111) strained case [8] and the number of electron pockets 
changes from six to three. When the band with double 
minima changes into that with a single minimum as a 
function of one parameter, it can be analytically shown, 
using the same treatment as Landau's phenomenologi- 
cal theory for second phase transitions, that the curva- 
ture at the minimum becomes zero; that is, the effective 
mass along the merging direction diverges. Consider- 
ing the time reversal symmetry, the band with double 
minima near the BZ boundary along the reciprocal lat- 
tice vector is expressed to the lowest order of k — ko as 



e(fc,A) = a(A) + b -^(k - k ) 2 + ^(k 



c(A), 



A^Ao 



/co) 4 , where A 

is the parameter used to change the band structure, k is 
the position in the BZ along the vector, and ko is at the 
BZ boundary. The k points with the extrema are given 

by k = ko , fco ± \/^xj • Assuming that, when A — > Ao, the 
double minima merge into one, the following condition is 
obtained: 



= 0, that is, b(X)\ x ^ Xo = 0. This 

means that the second derivatives of the band energy 
vanish and the effective mass diverges at Ao- 

The mass divergence occurs when the band structure 
with the double minima changes into that with the single 
minimum. This band structure change is derived from 
the strain and confinement. 

About the strain for the bulk, it is reported by Ma et 
al. that the degeneracy of the lowest conduction bands at 
the X point splits due to the (110) and (111) strains but 
does not do so due to the (001) strain as determined using 
a semiempirical tight-binding method and group theory 
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FIG. 5: Effective mass at A- valley and minimum k-position of 
the (110) strained Si thin film and bulk. Longitudinal effective 
masses of Si slab with 65 Si layers (~ 12nm) (closed square) 
and of Si bulk (open square). The k-position where the low- 
est conduction band takes the local minimum along the Y-X 
line ([001] direction) for Si slab ("+") and for bulk ("x"). 
The scale of the k-position is represented by the fractional 
coordinates, whose unit is FX (half of the nearest neighbor 
reciprocal lattice vector along [001]). The Brillouin zone (BZ) 
of the slab is folded and the A-point at the BZ boundary of 
the bulk corresponds to the T-point in the slab. To allow it 
to be drawn in the same figure, the k-position for the slab 
is transformed by 1 — k/2, where k is the k-position for the 
slab. The bulk effective mass along the confinement direction 
([110]) at the k-positions of 0.84 (closed circle) and 0.95 (open 
circle) . 



[2lL [22t . In this study, it is shown that sufficiently strong 
(111) and (110) strains largely break the symmetry, split 
the degenerated bands at the BZ boundary, and make the 
lower band with the single minimum at the BZ boundary, 
which causes the effective mass anomalies for the thin 
films and the bulk. 

The effect of the confinement on the effective mass 
anomalies is qualitatively explained as follows. Accord- 
ing to the effective mass approximation (EM A), the effect 
of confinement on the band structure is described by the 

h 2 2 1 

increase of zeropoint energy 2 m f > wnere m conf and 
L are the effective mass along the confinement direction 
and the thickness, respectively. We extend this EM A ex- 
pression using the effective mass at the energy minima to 
the whole band near the BZ boundary. In Fig. [5j the ef- 
fective mass along the confinement direction [110] at the 
k-points of 0.84 and 0.95 in the fractional coordinates is 
plotted using closed and open circles, respectively. The 
value 0.84 corresponds to the local minima without strain 
and 0.95 to the neighborhood of the BZ boundary. In the 
tensile strain region(> 1), the mass at 0.95 is larger than 
that at 0.84. Therefore, as the confinement thickness de- 
creases, the energy level at the 0.84 increases compared 
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with that at the BZ boundary and the structure with 
the single minimum is preferable to that with the dou- 
ble minima, which well explains the behavior of effective 
mass anomalies in Fig. [2j In the compressive region, the 
same discussion is valid. A similar discussion is found to 
hold for the (111) confinement case. 

For the experimental observation, it is in the (110) 
confinement case that the effective mass anomalies have 
the largest effect on the electron transport, because, in 
the (111) confinement, the anomalies are averaged out 
by the equivalent valleys along three different directions. 
For the (110) confinement, there are two kinds of valleys: 
A and A', in which the latter does not reveal the anoma- 
lies. The parameter region suitable for the observation 
corresponds to the anomalies where the A-valley is the 
conduction bottom. In the (110) confinement, there are 
two such regions: that of the thickness range from the 
bulk to 6 nm and around the 0.98 strain, and that of 
the thickness range less than 2 nm and the strain range 
from 0.99 to 1.01. It is experimentally reported that the 
mobility for channel thicknesses less than 5 nm decreases 
in the UTBD FET with the (110)-confined[6]. This mo- 
bility degradation may be due to the effective mass en- 
hancement as well as the surface roughness scattering. 
When the conduction current is large, the anomaly points 
where the A-valley is energetically higher than the A'- 
valley may affect the transport, because the energy differ- 
ence between the A- and A'-valleys is small (<50 meV). 
For applications that involve, combining piezoelectric de- 
vices, the effective mass may be controlled by voltage. 

Note that, considering the mechanism of symmetry 
breaking and confinement, the effective mass anomalies 
are expected to be found for materials with a similar band 
structure to Si such as diamond (C), SiC, AlSb, and GaP. 

In summary, the effective mass of strained Si thin films 
and crystals was systematically investigated using first- 
principles calculations. In (111) and (110) strained Si, 
effective mass anomalies were found for both the thin 
films and the crystals. The confinement and the symme- 
try breaking due to the strain changed the band struc- 
ture from double minima to a single minimum, which 
caused the anomalies. In such a case, it was analytically 
shown that the effective mass diverges. The dependences 
of the confinement on the anomalies were qualitatively 
explained by extending the effective mass approximation 
to the whole band near the BZ boundary. The experi- 
mentally observable regions of the strain and thickness 
were obtained for the (110) confinement. 
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